Hall Coefficient of Dirac Fermions in Graphene under Charged Impurity Scatterings 
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With a conserving formalism within the self-consistent Born approximation, we study the Hall con- 
ductivity of Dirac fermions in graphene under charged impurity scatterings. The calculated inverse 
Hall coefficient is compared with the experimental data. It is shown that the present calculations 
for the Hall coefficient and the electric conductivity are in good agreement with the experimental 
measurements. 
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Since the experiments of graphene were realized [j], 0, 
3], much effort has been devoted to studying the trans- 
port properties of the Dirac fermions. Many theoretical 
studies are based on the model of zero-range scatters in 
graphene @, S, 0, B H E3, E3- However, it has been 
found that the charged impurities with screened Coulomb 
potentials [l2|, [H, 13| are responsible for the observed 
carrier density dependence of the electric conductivity of 
graphene (2|. Although the electric conductivity experi- 
ment has been successfully explained, so far there existed 
no satisfactory theories to fit the density dependence of 
inverse Hall coefficient as measured by another experi- 
ment in Ref. 0. 

In this work, with a conserving formalism within the 
self-consistent Born approximation (SCBA), the Hall 
conductivity is calculated by using the diagrams gener- 
ated from the current-current correlation function. Wo 
show that the experimentally measured electric conduc- 
tivity and the inverse Hall coefficient can both be suc- 
cessfully explained in terms of the carrier scatterings off 
charged impurities. 

At low carrier concentration, the low energy excita- 
tions of electrons in graphene can be viewed as massless 
Dirac fermions 0- 
cent experiments 



17| as being confirmed by re- 
. Using the Pauli matrices it's and 
t's to coordinate the electrons in the two sublattices (o 
and b) of the honeycomb lattice and two valleys (I and 
2) in the first Brillouin zone, respectively, and suppress- 
ing the spin indices for briefness, the Hamiltonian of the 
system is given by 
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) is the fermion operator, 



the momentum k is measured from the center of each 
valley, v (~ 5.86 eVA) is the velocity of electrons, V is 
the volume of system, and Vi{q) is the Thomas-Fermi- 
type charged impurity potential Here, we neglect 
the intervalley scatterings in Vi{q) for two reasons. First, 
for low electron doping, the intervalley scatterings are 
negligible small than the intravalley ones. Second, by 
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FIG. 1: (color online) (a) Self-consistent Born approxima- 
tion for the self-energy, (b) Current vertex with impurity 
insertions, (c) Diagrams for calculating the Hall conductivity 
a xy . The vertex j is associated with the vector potential Aj. 
k — k ± q/2. = io n + Q, m with Q m and uj„ respectively 
the bosonic and fermionic Matsubara frequencies. 



doing so, our formulation of the problem given below 
will be much simplified. 

Under the SCBA [see Fig. 1(a)] [H, 0, the Green 
function G{k,uj) = [to + fi — vk ■ <jt z — E(/c,w)] _1 = 
go(k,ui) + g c (k,Lu)k ■ <7T Z and the self-energy ta) of 
the single particles are determined by coupled integral 
equations [14| . Correspondingly, the current vertex cor- 
rection is the ladder diagrams in Fig. I (b). The current 
vertex vT x {k,u 1,012) is expanded as 

3 

T x {k,w u u 2 ) =^£ j yj{k,w 1 ,u) 2 )A x j {k) (2) 
j=o 

where ^4§(fc) = r z a x , Af (fe) = a x a ■ k, Af (k) = 3 ■ ka x , 
= t z 3 ■ ka x a ■ k, and j/jffc, u>i, UJ2) are determined 
by four-coupled integral equations [14j. 

We start to calculate the Hall coefficient of the Dirac- 
fcrmion system in graphene. Consider that the system 
is acted with an weak external magnetic field B per- 
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pendicular to the graphene plane. The vector poten- 
tial A(q) is related to B via B = iq x A(q) where q 
is a wave vector. The Hall conductivity a xy is defined 
as the ratio between the electric current density along 
x direction and the electric field applied in y direction. 



By the standard linear response theory, a xy is obtained 
as a xy = — limn-,o W^jfJl + z0 + )/f2 where <f> xy (SY) is 
the current correlation function. For the Matsubara fre- 
quency f2 m = 2-KmT with m as a integer, <f> xy (i£l m ) is 
given by [201 ] 



* a »(*n m ) = --V / dr / dT / e* n » r (T T J r a! (flf,T)J 1 ,(0 ) 0)J i (-«,T y ))i4 i ( 9 f) (3) 
^ j Jo Jo 

where T T is the r ordering operator, J x (q,r) is the current operator, /3 = T _1 , (• • • ) means the statistic average, 
and the use of units in which h = c = ks = 1 has been made. Within SCBA, the function <& xy (iQ, m ) is calculated 
according to Fig. 1(c). Writing it explicitly, we have 

$ xy (iQ m ) = 22 Tr {[T x (k~, k + , iu n , iu+)G{k + ,iu+)T y (k + ,k + ,iuj+ ,iuj n ) + T y (k~ , k~ ,iu> n , iuj~) 

knj 

xG(k~,iu~)T x (k~, k + iiu n )]G(k + ,iu n )Y j(k + , k~ ,iuj n ,iu; n )G(k~ ,iuj n )}Aj(q) (4) 

where the factor 2 stems from the spin degeneracy, k = k ± q/2, and u„ — w„ ± f2 m . The vertex given by Eq. ([2|) 
corresponds to T x (k, k, u>i, 0J2) = ^x(k,0Ji,u>2)- r p (fc~ , k + , oji, u/2) satisfies the 4x4 matrix equation 

r M (fc~, k + , wi,cj 2 ) = t 3 ct m + -i^riiUo(fc - fci)G(fcr, Wi)r M (fcf , fc+, wi, w 2 )G(fcj f , w 2 ). (5) 

fci 
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r /J (fc + , fc~,a;i,a-> 2 ) is obtained by exchanging - and + in 
Eq. ((5|). The analytical continuation iil m — > i7 + i0 + 
for <& sv (iri m ) can be manipulated according to the text 
book [2l|. For the weak magnetic field B, we expand the 



right hand side of Eq. ^ to the linear terms of q (giving 
rise to the linear terms of B). To do so, we need to note 
the facts listed below. (1) T^(k~ , k + , uji, uj 2 ) is expanded 
as r M (/c,wi,w 2 ) +7 M (fc, g,wi,w 2 ) with 



7^(^,9,^1,^3) = —^2n i vl(k-k')G(k',uJi)j fl (k\q,uii,uj2)G(k' ) uJ2) 
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n ^ k - fc ')[VG(fc', Wl )r p (fc', wi, u 2 )G(k\ lo 2 ) - G(k', wi)r„(fc', wi, w 2 )VG(fc', wa)] • <f/2 

fe' 

(6) 

where V means the gradient with respect to k'. (2) From the identity 

V £ + aF )w ° 2(fc ~ *')]^[G(A:, w 07m(*> 9, wx, W2)G(fc, u>2)G(V, LJ 2 )?u(k', U2,ui)G(k', ujx)] = 0, 

fefe' ^ j 

rewriting the left hand side by performing the integral by part and using the equations for T v (k, u>i, w 2 ) and 
7 M (/c,g,wi,w 2 ), we obtain 

^2 Tr {^(^' 0> w i> w 2)[^-G(/c, w 2 )r i/ (fc, w 2 , wi)G(fc, wi) + G(fc, w 2 )r„(fe, uj 2 , u 1 )—G(k, wi)]} 

— * o 

= -| • V Tr{[VG(fc, wi)r^(fc, Wlj w 2 )G(fc, w 2 ) - G(fc, wijr^fc, Wl> w 2 )VG(fe, W2 )]— r„(fc, w 2 , wi)}. (7) 

fe i 
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(3) Using Ward identity, vG(k, oj)T a (k, u, oj)G(k, oS) — dG(k,ui) /dk a , we have 



vG(k+ ,u;)T a (k+ ,k~ ,Lo,io)G{k- ,u) = G(k, u) + \ ■ ^ ^{k)[a 3 {k, + ia z bj{k, u)4>], (8) 



3 

where <f> is the angle of k and <j> is the unit vector in <f> direction. This expansion can be further simplified. If 
we consider the left hand side as a functional of the matrices AJ(k), then its transpose should be a functional 

of matrices Af(k) defined as Ag(k) = T z a a , A?(k) = a^a 1 ■ k, A%(k) = a 1 ■ far*, and A$(k) = t z o 1 ■ ka^a* ■ k 
where a means the transpose of the Pauli matrix. By comparing the transpose of Eq. {SJ and the expansion of 
vG t (k~,uj)T t a (k + ,k~,uj,uj)G t (k + ,uj) in terms of A"(k), we find a\ = —a-} = a, b\ = 62 = b, and all other a's and 6's 
vanish. Using <j a a ■ k — a ■ ka a = —2ia z a ■ 4>, and a a u ■ k + a ■ ka a = 2a ■ k, we obtain 

Q .... 
vG(k + , uj)T a (k + , k ,aj,u))G(k ,u>) — — — G(k, ui) + ia z a ■ [b(k, u:)k(f> — a(k, uS)4>k] ■ q. (9) 

oka 

This result can also be obtained by solving Eq. |6]). Since the final result depends on a(k,uj) + b(k,u) = c(k,u>), we 
here present only the equations for determining c(k,uj) = z(k,ui) + [<7o(fc, u>) — g^(k,u))]X(k,u>): 

z(k,uj) = [g' (k,uj)g c (k,u;) - g (k,uj)g' c (k,uj)}[yo(k,uj,uj) - y 3 (k,uj,uj)} 

- 9c{k, uj){g (k, u)[y (k, ui, cj) + y 3 (k, u, w)] + 2g c (k, uj)yi(k, w, w)}/fc, (10) 

X(k,u) = i^ni^(* - fc'JWfc'.w) + fagtfc'.w) - ^(f.wp^.w)}, (11) 

k' 

where g'(k,u>) — dg(k,ui)/dk. Using above results, we obtain for the Hall conductivity 

Bv 2 e 3 ^ f°° dtu 8F{uj) ^ r r -+, n+( dT V dG- dT+- dG~ 



>xy 



V j^oa 2it duo x dk x dk y dk y dk x 

+ C^-w^- 2 '*- )G+r » + ^ i) (i2) 



where T x + = T x (k,uj ,w+), T+ = T y (k,uj + ,Lu ), 
and G ± = G(k,u> ± ). Some terms such as 
T-+{dG+/dk x T+-dG~/dk y - dG+ / dk y T+- dG~ / dk x ) 
and those containing the same frequency-w + arguments 
are not included because they happen to be zero under 
the operation ImTr in the right hand side of Eq. (Tl3>|) . 

The numerical result for the Hall conductivity a xy as 
a function of carrier concentration S (doped carrier per 
carbon atom) is shown in Fig. 2. Here the impurity den- 
sity Ui — 1.15 x 10~ 3 a~ 2 (with a as the lattice constant 
of graphene) is chosen as the same as in our previous 
work 14]. In the normalization factor 1/aB given in 
Fig. 2, the magnetic field B is a constant and a is the 
electric conductivity (see Fig. 3). In the limit S — ► 0, the 
normalization denominator is a constant because of the 
minimum conductivity. While at large carrier concentra- 
tion, it is linear in S because of a oc 5. Within a very 
narrow range of S around 0, a xy varies dramatically and 
vanishes at 6 = 0. Beyond this regime, the saturation 
behavior of the curve implies that cr xy i pc 8. The present 
result differs quantitatively from Ref. ly using a different 
approach. Also, it is qualitative different from Ref. [22j 
in which a constant scattering rate was phenomenologi- 



cally introduced and the current correlation was treated 
without vertex correction. 

For calculating the Hall coefficient, we need the result 
of the electric conductivity a obtained with the same 
scattering parameters as for a xy . The result for a as a 
function of 6 is shown in Fig. 3 and compared with ex- 
perimental data |2|. The minimum conductivity by the 
present calculation is about 3.5 e 2 /h close to the data 4 
e 2 //i, which is a consequence of the coherence between 
the upper and lower band states. Our numerical calcu- 
lation shows that the contribution from intervalley scat- 
terings is negligible small. 

In Fig. 4, we exhibit the theoretical result (solid line) 
for the inverse Hall coefficient defined as R^ 1 — Ba 2 /a xy 
and compare it with the experimental data (sym- 
bols). Clearly, the present calculation fits the exper- 
imental measurement of the inverse Hall coefficient as 
well as the electric conductivity very well. With compar- 
ing to the classical prediction R^ 1 = —nec (with n as 
the doped electron density) , both the present calculation 
and the experiment data for diverge at 8 = 0. The 
divergence of R^ 1 stems from the vanishing of a xy at 
5 = while the conductivity a remains finite. The clas- 
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FIG. 2: (color online) Hall conductivity a xy normalized by 
aB as function of the electron doping concentration 5. a is 
the electric conductivity, and B is the magnetic field in unit 
of Tesl; 
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FIG. 4: (color online) The inverse Hall coefficient (in 
unit of 10 _3 Tesla/Ohm) as function of the electron doping 
concentration 8. The present calculation (solid line) is com- 
pared with the experimental data in Ref. |2| (symbols). The 
dashed line represents the classic theory i? -1 = — nec. 



at zero carrier concentration stems from the vanishing of 
the Hall conductivity and meanwhile the minimum re- 
mained in the electric conductivity. The present results 
for the inverse Hall coefficient and the electric conduc- 
tivity are in very good agreement with the experimental 
measurements. 

This work was supported by a grant from the Robert 
A. Welch Foundation under No. E-1146, the TCSUH, the 
National Basic Research 973 Program of China under 
grant No. 2005CB623602, and NSFC under grant No. 
10774171 and No. 10834011. 



FIG. 3: (color online) Electric conductivity a as function of 
the electron doping concentration 5. The present calculation 
(solid line) is compared with the experimental data in Ref. [2J 
(symbols). 



sical theory is based on the concept of the drift velocity. 
At the zero carrier concentration if the conductivity re- 
mains finite, the drift velocity has to become infinitively 
large, which implies an infinitely large Lorentz force act- 
ing on an electron. Therefore, the classical theory is not 
applicable near 5 = 0. On the other hand, at large car- 
rier concentration, the present calculation reproduces the 
classical theory, and both of them are in agreement with 
the experimental result. 

In summary, on the basis of self-consistent Born ap- 
proximation, we have calculated the Hall coefficient of the 
Dirac fermions under the charged impurity scatterings in 
graphene. The anomalous in the inverse Hall coefficient 
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